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'sj- ' Abstract 

I We study spin glasses on random lattices with finite connectivity. In the infinite 

' connectivity limit they reduce to the Sherrington Kirkpatrick model. In this paper we in- 

. vestigate the expansion around the high connectivity limit. Within the replica symmetry 

breaking scheme at two steps, we compute the free energy at the first order in the expan- 
sion in inverse powers of the average connectivity (z), both for the fixed connectivity and 
for the fluctuating connectivity random lattices. It is well known that the coefficient of 
I the 1/z correction for the free energy is divergent at low temperatures if computed in the 

Q ' one step approximation. We find that this annoying divergence becomes much smaller 

O ' if computed in the framework of the more accurate two steps breaking. Comparing the 

temperature dependance of the coefficients of this divergence in the replica symmetric, 
KM \ one step and two steps replica symmetry breaking, we conclude that this divergence is 

^ • an artefact due to the use of a finite number of steps of replica symmetry breaking. The 



1/z expansion is well defined also in the zero temperature limit. 
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1 Introduction 



Many studies have been devoted to finding analytic solutions of more realistic models than the 
Sherrington-Kirkpatrick one. The diluted spin glass models belong to this class, and they are 
characterized by a finite coordination number; these models are also interesting because they 
are connected with different optimization problems 

In the present work we consider lattices where each site is connected with a finite number 
of randomly chosen sites; we study both the cases where the connectivity is fixed and where 
the connectivity is a Poissonian variable with given mean value. The spin interaction is only 
among nearest neighbour pairs. 

The random structure of these lattices allows us to neglect the probability of closed paths of 
finite length: this probability becomes indeed zero in the thermodynamic limit: the correlations 
among the neighbours of a given spin can be neglected. We are therefore dealing with mean 
field models, although the difficulties due to the finite connectivity don't allow us to solve them 
exactly. 

Whereas in the SK model only the overlap between two replicas occurs as order parameter 
(the order parameter is a function in the infinite step replica symmetry breaking solution and 
a pure number when the replica symmetry is exact), in the finite connectivity models the order 
parameter becomes a function of the overlaps of any number of replicas and then it becomes 
a function of an infinite number of variables when the symmetry is totally broken; as a result 
it is extremely difficult to find the exact free energy 0, |^, In other words the probability 
distribution of the effective cavity fields is Gaussian in the SK model as a consequence of 
the central limit theorem, so it can be characterized by its variance. When the number of 
neighbours z is finite, this distribution is no more Gaussian and all the moments are relevant 
and this leads to the presence of an infinite numbers of order parameters (also in the replica 
symmetric situation). 

Perturbative solutions have been investigated both near the critical temperature [§] and 
near the infinite connectivity point (SK model) 0, 0. Recently it has been proposed 
proposed a general non perturbative solution developing the Bethe-Peierls cavity method to 
an approximation that is equivalent to a one step replica symmetry breaking level 0. 

The present work addresses to the study of the large connectivity expansion: we compute 
the first order of the expansion in the inverse power of the connectivity [z) for the free energy. 
The 1/z expansion has been studied for the fixed connectivity model by Goldschmidt and De 
Dominicis, at the first step of replica symmetry breaking |]^; they found results that exhibit 
a low temperature divergence for the first order correction in 1/z to the free energy density. 
The 1 / y/z expansion, that they computed a T = and to the second step of replica symmetry 
breaking, has finite corrections both at the first and the second order [Q. The Ij \fz coefficient 
becomes yet smaller by a factor ten and by a factor three when one goes from the replica 
symmetric solution to the IRSB one and from this to the 2RSB one respectively. 

In this paper we have computed the coefficient of the Xjz expansion up to the second step 
of replica symmetry breaking. For the replica symmetric and the IRSB solutions our results 
agree with those found by Goldschmidt and De Dominicis. These results suggest that the 
pathological behaviour (i.e. the low temperature divergence) is a consequence of the fact that 
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one stops the computation at a finite step of tlie iterative process for breaking the rephca 
symmetry. Our results indicate that the 1/z expansion is well defined and can be used also 
in the zero temperature limit. We notice that a well defined 1/z expansion is possible for a 
model with continuos varying coupling only if (irrespective of the sign) the zero temperature 
entropy is zero in the limit z ^ oo. Indeed it is easy to prove using the approach of that in 
the mean field approximation for finite z the zero temperature entropy is identically zero, so 
that the two limits z ^ oo and T — > could not be exchanged in an hypothetical model with 
continuos coupling if the zero temperature entropy were different from zero at z = oo. 

The paper is organized as follows: in the second section we present the two models we 
study and the high connectivity expansion is obtained. We use a simple method to evaluate 
sums over multiple replicas overlaps. In the third section we show how to perform sums over 
the replica indices in a simple way and in section four we illustrate our numerical results for 
the two steps replica symmetry breaking and we compare them with the known ones at the 
first step of replica symmetry breaking. Finally we present our conclusions. The appendix is 
devoted to a consistency check for the form of the free energy we use. 



2 The large connectivity expansion. 

By definition a Bethe lattice is a lattice where the Bethe-Peierls approximation is exact; this 
is equivalent to saying that there are no finite size loops. 

In the random lattices we study the typical length of a loop is proportional to logA^: in 
the infinite volume limit it is therefore a Bethe lattice. This is locally equivalent to a tree- 
like structure, nevertheless by defining the Bethe lattice as a random lattice one bypasses the 
problem of fixing the boundary conditions to introduce frustration (this is provided by the 
loops of size ~ In A^) . 



2.1 Random lattice with fixed connectivity 

We therefore follow a variational formulation using in the framework of the replica approach 
the same scheme |^. We define a functional and we show that the free energy is obtained as 
the stationary point with respect to an order parameter that will be defined. The free energy 
functional is 0]: 

n/3/„((7„) = zln(Tr{.4^7r'({aJ))+ (1) 




+ 00 



rfJP(J)Tr|.4Tr{.4(7:({aJ)^7:({r4)exp[/3J5:a, 



a=l 



where Trs^^j^^ is the sum over the 2" configurations of the variables cxa with a = l,---,?T,,2 + l 
is the lattice connectivity, 5'n({o"a}) is a function of the n variables aa and it plays the role 
of the order parameter. Our goal is to make stationary the functional /„ with respect to the 
variation of 5'n({o'a})- We have then to find the solution of the equation: 

f = 0. (2) 

5gn 
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that gives for the order parameter the equation: 



gniWa}) = C / dJP{J) Y: exp [Y.^J^aTa h?^({r„}) (3) 

{r4 \a=l J 



^ r + OO .] jT,f T\rn.. rn,. r _ T N r_ T N to T v^ra _ _ 1 • 



with 

dJPiJ)Tr{^^yTr{,^yg^,{{a,})g^^{M) exp [/^JELi ^a^. 

We can notice that the functional defined by (|l]) is independent of the Qn normahzation; we 
can then use a convenient one, provided one changes the constant C (in (|^)) in CS-^~^^ when 
changing Qn in Qnd. The correctness of this functional (0) has been proved by De Dominicis et 
al. ; a simple way to get this result is reported for completeness in the appendix. 

In order to write the order parameter gn{{<^a}) in a more explicit form (where the multiple 
overlaps appear) we generahze the identity 

exp (pjaaCTb) = cosh (/5J)(1 + aaCTfo tanh {(3 J)) (5) 

to 



exp ( /? J aaCiA = cosh" (/? J) J2 ( tanh^ (/? J) ^ - 

\ a=l / ■r=0 \ ai<-- <ar 



(6) 

where the last sum is over all possible sets of r replicas, counting once any permutation. 

2.1.1 Interaction with a bimodal distribution. 

We first study the following distribution for the J: 

PiJ) = l[S{J + Jo) + SiJ-Jo)] . (7) 

Equation (^ is formally identical after averaging on the J, providing one sums only over the 
even r and writes Jq instead of J. 
If we define the overlaps 

we can write the eq. (|) as: 

9n{Wa}) = cosh" {(3 J) J2 ftanh^ {f3J) ■ ■ ■ aa.ga,..aJ • (9) 

r=0 \ ai< - <ar / 

We can now implement the ^ expansion if we scale the couplings as usual: 

J =4= (10) 
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and set J = 1. Performing the expansion, after some computations we obtain at the first order: 



with: 



and 



/ = /o + -/i + 0(4), (11) 
/5/o = + f E ^if - - HTr exp(/5 ^ ^i^^a^f )] (12) 

/^/l = -T + ^-^^^-^^22<lab E 5a6cd (13) 



+ Wa6 Ibc lea ) ^ 2^ [lab led + lae %d + ybc ) labcd ■ 

a<b<e a<b<e<d 



As it should be (we are expanding around z — +00), /o is the SK free energy. In these 
expressions we have also expanded the overlaps in powers of 1/z: 

- «(") J- ^n(^) 4- 
Qab - Qab + -Qab "I ' 

qabed = Qabed + l(labed + --- (14) 
Z 

and we have used the identities: 

(0) Tr^exp[p^J2r<sQi^^'^r(^s]craCrb ^ ^ 
^ab = -— r— (0) —=<(Taa,>Q a^b 



(0) Tr^ exp /j^ ^ qf}(JrCrs](ya(Tb(Te'7d _^ ^ / ?. / / J ^ ^^ 

Tabed = T— ; m : =< (7a(Jb(7eCrd >Q a ^ h ^ C ^ d , (15) 

Tr^ exp [(3^ Y.r<s Qrs'arasl 

where < ■ >q is the average on the single site Sherrington-Kirkpatrick Hamiltonian. We 
notice that / is no longer stationary with respect to the order parameter q because we have 
already used the stationary equations to simplify the result. 

2.1.2 Interaction with a Gaussian distribution. 

If we use a Gaussian distribution with the same mean (J = 0) and variance {Jq/ z) of the 
previously studied bimodal distribution, one finds that at this order the only relevant difference 
is in the fourth moment of the interaction (SJq/z^ for the Gaussian and Jq/z"^ for the bimodal 
one), (it is crucial that at this order we expand the order parameter only up to the second 
order in z, i.e. up to the fourth order in J). Performing the same calculations as before, we 
arrive to the final form for the free energy first order correction: 

a<b a<b<c<d 
a<b<e a<b<e<d 
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As expected, the /o does not change because it does not contain terms and the SK model 
is indeed independent from the particular distribution one uses if we fix the mean and the 
variance of the couphngs. 



2.2 Random lattice with fluctuating connectivity. 

In an other interesting model the connectivity is a Poissonian variable with mean value z. 
We take into consideration the large z expansion, where the interactions probability distribution 
can be written in the form: 

P( J,,) = (1 - ^)5{J,k) + ^P{J^k) Vz, k , (17) 

where P{Jik) is a distribution to be defined. 

In principle we could write an expression similar to eq. (|l]) for the free energy, however it 
is simpler to proceed in a direct way. The n replicas partition function is: 

/ + 00 
P{J^k)dJ,kTr^ exp(/5 J,fe ^ «) (18) 

(z z r+°° ~ \ 

2.2.1 The expression of the free energy 

In the case where: 

P{Jik) = l[5iJ^k - Jo) + 5iJ,k + Jo)] Vz, k (19) 



we obtain: 



= l[TrJl + ^JcosHPJoY.<^t)-^]] 

i<k V a ) 



= Tr. exp{ ^^[cosh(/3 Jo E^^fc) -I]}- (20) 

i<k a 

Let us rescale: 

Jo - 4 (21) 

'z 



and write Jq = \. We can than perform the Xjz expansion up to the first order for the free 
energy: 

^ = ^-^exp{^5:[g(5:ax)V^(5:ax)']}= (22) 

i<k 0- 0- 

= Tr^ exp{l Ety (E ^t^>yk) + ^( E ^t^>'^l<yt<^l<^f<yt)]} ■ 

i<k a.b a,b,c,d 
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After converting the summations over replicas indices into distinct indices summations (using 
{(t"')^ = 1), introducing the Gaussian integrals and solving with the saddle point method, we 
find at the first order: /?/ = /3/o + ^Pfi, where /o is the SK free energy and /i has the form: 

a<b a<b<c<a 

In a similar way, when P{J) is a Gaussian distribution, we find the same result as before 
with the difference that the first order free energy /i is multiplied by a factor three. 
If we put together the previous formulae we find that 

Pfi = Afr'"" + Kf( (24) 

where A = 1 for the model with fixed number of neighbours, A = for the model with 
fiuctuating number of neighbours and K is the kurtosis of the distribution of couplings J. The 
quantity // is given by eq. (EBI) while y"*^*^'* by the formulae and O is found to be: 



= -i^-C(i-'^')E^lf (25) 



4 2n ^, 

a<b 



/ . \Hab Hbc ^ca ) + ^ \Hab Hcd + Hac Hbd + Had Hbc JHabcd 



^ a<b<c ^ a<b<c<d 



3 Evaluation of the sums over replica's indices. 

If the replica is broken at two steps (using the usual conventions) we can write: 

-Y.1ab = - l)ql + {mi - m2)ql + miql , (26) 

where the sum is on the indices a ^ h. Similar expressions can be written at higher orders in 
the replica symmetry breaking and in the continuum limit one obtains: 

hm-^g^ l\2^^y (27) 

n-.o Jo 

The formulae for the case where the replica symmetry is broken at two steps can be obtained 
by using: 

q{x) = qo < X < nil (28) 

q{x) = qi mi < X < m2 (29) 

g(x) = q2 m2 < X <1 (30) 



The direct computation of 



1^0n(n-l)(n-2)(n-3)5/'-^ ^''^ 
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Figure 1: Four replicas diagrams. The number below the trees is the degeneration due to the 
replicas indices permutations. 



is more involved and if it is not properly done it can become a nightmare. 

We can simplify it if we remark that the four replicas overlap is a function of all the possible 
x's among the four replicas: 



Qabcd q{xi2, xis, Xi4, X23, X2A-, 3:34), 



(32) 



with < Xij < 1 Vz, j; however, due to the ultrametric structure of the states at most three of 



the two replicas overlaps can be distinct. Mezard and Yedidia |T21 has shown that in order to 
compute this kind of sums over replicas is suitable to consider the five possible ways in which 
four replicas can be organized (fig. |l|); we have to associate a variable Xi to each vertex and a 
factor xl~'^{s — 2)! when s lines converge to it. Using this rule, providing to take into account 
the number of different permutations of replicas indices that produce the same configuration, 
we finally obtain: 



lim TT TT r 

n^o n{n -l){n- 2){n - 3) 



/■O px^ rX3 

y2llbcd = ^ dx2 I dxiq{xi,X2,XzY + 



X2 



dxiq{xi,X2,XzY + I dxi2x\q{x\) + 



rO r-xz 

+ 12 J dx^ J dx2 

/O rX2 rO rX2 

dx2 J Xidxiq{xi,X2Y + 6 j X2dx2 j dxiq{xi,X2) 

At the second step of replica symmetry breaking we obtain: 
f (llbcd= im2-l){m2-2){m2-3)ql^ 



(33) 



a<b<c<d 



3(1 - ma) (mi - m2)q2,223b - 3mi(l - ma) ^2222 



bd 



8 




3 1 

Figure 2: Three replicas diagrams. 

- 6(1 - m2)(m2 - mi)(2m2 - ^1)^222156 - 6(1 - m2)(m2 - mi)miql^2^bd 

- 12(1 - m2)mlq^^2o + ("^1 ~ "^2)("^i - 2m2)(mi - 3m2)ql^ 

- 3(m2 - mi)^mig2i2i - 12(m2 - mi)mlq2^2o - ^^illo 

+ 4(m2 - l)(m2 - 2)(mi - m2)ql^i^ - 4mi(m2 - l)(m2 - 2)^331^ 

- 4mi(mi - m2)(mi - 2m2)g3ji„ - 12mi(l - ni2){m2 - "^i)g22iiio (34) 



^ X! QabcdiQabflcd + lacQbd + QadQbc) — 
a<b<c<d 

3(m2 - l)(m2 - 2)(m2 - 3)^42^2 + 3(1 - m2f{mi - m2)q2222sb{ql + '^qD 

- 4mi(l - m2fq2222bd{ql + 2go) - 6(1 - m2){m2 - mi)(2m2 - mi)q2.,2^sb{q2qi + 2gi) 

- 6(1 - 7712) (m2 - mi)mig222iM(g2gi + 2go) - 12(1 - m2)mlq2.,2o{q2qo + 2go) 
+ 3(mi - m2)(mi - 2m2)(mi - 3m2)q4,ql - 3(m2 - mi)^mig2i2i (?? + 2go) 

- 12(m2 - mi)mlq2^2o{qiqo + 2go) - 18"^i94o?o 

+ 12(m2 - l)(m2 - 2)(mi - m2)g32ii'?2gi - 12mi(m2 - l)(m2 - 2)g32iog2go 

- 12mi(mi - 7712) (mi - 2m2)g3iio9i'?o 

- 12mi(l - m2)(m2 - mi)g22iiio(92go + 2gigo) , (35) 

where in the notation g^,, -- the quantity A is the number of rephcas and a indicates the block 
to which they belong (we refer to the matrix Qab at the second step of the ultrametric Ansatz); 
when two possibilities can occur, we write sb when the four replicas are in the same block of 
first replica symmetry breaking and bd in the other case (i.e. q222isb means that two replicas 
belong to the same second replica symmetry breaking block and all the four to the same block 
of first replica symmetry breaking; whereas q222ibd means that two replicas belong again to the 
same second replica symmetry breaking block, the other two to the same first replica symmetry 
breaking block but the overlap between the first two and the second two is the minimum one). 
For the sum on three replicas indices we have (see fig.(0)): 

! E (^i^M) = ("^2-l)(m2-2)g3 

a<b<c 

+ 3(1 - 7712) (m2 - nii)q2ql + 3(1 - m2)mig2'?o 

+ {mi - m2){mi - 2m2)ql + ^mi{m2 - mi)qiql + 2'm\ql . (36) 
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Substituting these expressions into 0, [T^, |2^, we obtain the exphcit expression for the free 
energies. We can now find the numerical values of g2, Qi, Qo, ""^2, f^i maximizing /o (the 
SK Hamiltonian is stationary with respect to q{x)) and then we can use these values in the 
expressions of the four replicas overlaps. 

To obtain the expressions at one level of RSB we can put q2 = qi and identify the four 
rephcas overlaps in this way: q^^ = q^^ = g2222sb = Q222isb = feu; ^Saio = Qsiio = feiiio! 

q2222bd = Q'2i2i = q222ibd', ^20 = 5'2i2o (sGC [0]). 



4 The solution of the equation 

To evaluate the value of the free energy we have firstly to solve for the m and q parameters 
of the infinite connectivity limit and to compute the parameters in equations (pT]), (p^D, (pS]), 
56]). At this end we have to compute integrals like the following: 



Q2222sb 



dz , 



Num 
Den 



Num 



(i?/exp( 



dx exp ( 



2(gi - go) 




X 



2(g2 - qi) 



dx exp( 



X 



2(g2 - qi) 



cosh'"2(/3(2 + ?/ + x)) 



tanh^(/?(2 + y + x)) cosh"'"(/5(2 + y + x)) 



Den 



(i?/exp( 



2(gi - go) 




dx exp(- 



2(g2 - gi) 



cosh'^^P^z + y + x))) (37) 



To evaluate these expressions it is important to optimize the number of operations the 
computer has to do. In the numerical evaluation of the integrals (we are considering sums 
instead of integrals and we set x = a z where i is an integer) the internal integrals have in fact 
to be evaluated for every value of the variable of the external one. A repeated evaluation would 
take an enormous amount of time. A much faster method consists in evaluating beforehand 
the internal functions (i.e. cos\i{x + y + z)) for all values x + y + z = ai and in storing in a table 
the values in the integrals; in this way the computer has to perform a number of operations 
proportional to instead to A^'^ of the naive method (A^ is the number of spacings in which 
the integration domain is divided). 

The final results for the coefficient of the 1/z corrections to the free energy are shown in 
fig ^ in the case of ±1 interactions. We immediately see that the divergence of the correction 
to the free energy at T = fades away when we increase the order of the replica breaking and 
it is an artefact of using a starting point which is not correct (the correct one corresponds to 
infinite breaking of the replica symmetry). 

We evaluated the entropy doing the derivative of the free energy {S = —df/dT) using an 
high order expression for the finite difference derivative. The final results for the coefficient 
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Figure 3: The Xjz correction to the free energy as function of the temperature for the 
rephca symmetric case (*), one step rephca symmetry breaking and two steps rephca symmetry 
breaking for the model with J = ±1 for fixed connectivity (left) and for fluctuating connectivity 
(right), 



-0.26 
-0.27 
-0.28 
-0.29 
-0.3 
-0.31 



+ « 
+ 



2RSB 
IRSB 
RS 



-0.32 I ' ' ' ' 1 

0.02 0.04 0.06 0.08 0.1 

j3 



Figure 4: The 1/z correction to the free energy as function of for the replica symmetric 
case (*), one step replica symmetry breaking and two steps replica symmetry breaking for the 
model with J = ±1 for fixed connectivity 
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0.15 
0.1 

0.05 



2RSB 
IRSB 
h(x) 
hl(x) 



0.1 0.2 0.3 0.4 0.5 

T 



Figure 5: The 1/z correction to the entropy as function of the temperature for one step rephca 
symmetry breaking and two steps rephca symmetry breaking for the model with J = ±1 for 
fixed connectivity and the corresponding second order polynomial fits in the high temperature 
region. 

of the 1/z corrections to the free energy are shown in fig ^ in the case of ±1 interactions 
with fiuctuating connectivity. Also in this case we see that the divergence of the correction 
to the entropy near T = fades away when we increase the order of the replica breaking. 
The correction for the entropy are much stronger that those for the free energy. In order 
to evidentiate the effect of the spurious divergence at T = we show also a second order 
polynomial fit in the high temperature region, which dramatically fails at low temperature. 



5 Conclusions. 

We can see from the numerical data that the first order correction fi of the free energy in 
all the analyzed models taken into account presents a divergence at small temperatures. We 
fit in the interval T G [0.05, 0.5] a behaviour of the kind (see fig.(|^)): 



C 

f^(T) =D + AT^ + BT^ + — . 



(3^ 



We report in the following table the values of the C coefficient in the different cases considered. 
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C values 




fixed connectivity 


fiuct. connectivity 


±1 


Gaussian 


±1 


Gaussian 


RS 


0.035 


0.12 


0.35 


0.12 


IRSB 


0.003 


0.010 


0.003 


0.010 


2RSB 


0.001 


0.003 


0.001 


0.003 



This coefficient is three times smaller when going from the IRSB solution to the 2RSB 
one. Moreover, if we look at the results in we can see that the same ratios have been found 
for the 1/v^ coefficient; its values are in fact 0.1 for RS, 0.01 for IRSB and 0.0026 for 2RSB 
solutions. The same divergent factors occur in the fixed and in the fluctuating connectivity 
models and the entire divergence comes from the term proportional the forth moment of the 
distribution of the J. However, we can see from the figures that the divergence moves to 
smaller temperatures when the number of replica symmetry breaking steps increases. In the 
fixed connectivity model with bimodal distribution the divergence appears at T < 0.1. 

We notice that in all these models there is a correction to the low temperature behaviour 
of the S.K. model, where the free energy is proportional to T^. From the numerical data suggest 
that this correction goes to zero in the full replica symmetry breaking solution, at least in the 
fixed connectivity model, that is linear in over a large range of temperature already in the 
2RSB solution (fig.(|)). 

We ffited the curves in the temperature interval where the divergence doesn't yet occur, 
with the function: 

f^{T) = D + AT^ + BT'^ . (39) 

For completeness we report in the following table the values of the coefficients A and B in the 
different cases considered. 



B values 




fixed connectivity 


fiuct. connectivity 


±1 


Gaussian 


±1 


Gaussian 


RS 


-0.199 


-0.049 


-0.63 


-1.88 


IRSB 


-0.122 


-0.139 


-0.37 


-1.12 


2RSB 


-0.054 


-0.08 


-0.3 


-0.9 



A values 




fixed connectivity 


fiuct. connectivity 


±1 


Gaussian 


±1 


Gaussian 


RS 


-0.286 


-0.156 


0.19 


0.619 


IRSB 


-0.30 


-0.176 


0.165 


0.497 


2RSB 


-0.313 


-0.18 


0.16 


0.48 
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The fit we have done assumes that, apart form the divergence, the entropy extrapo- 
lates to zero at zero temperature^. In order to check the consistency of the results, we have 
extrapolated to zero temperature the numerical results for the entropy in the temperature 
interval where the divergence doesn't appear yet. We then find a behaviour that accords with 
the expected one: the zero temperature value s(0) is different from zero but decreases when 
going to higher steps of replica symmetry breaking, suggesting that it will reach the correct 
value Si(0) = in the infinite steps limit. 

We give the results in the following table. 



1st order correction to entropy at T = 




fixed connectivity 


fiuct. connectivity 


±1 


Gaussian 


±1 


Gaussian 


IRSB 


0.26 


0.48 


0.12 


0.36 


2RSB 


0.18 


0.32 


0.07 


0.22 



To conclude, we think that there are numerical evidences that confirm that the Xj z expan- 
sion is correct to study the random lattices with high connectivity. The Xjz expansion arise 
naturally if we compare the high connectivity limit in random lattices to the case of nearest 
neighbor interactions in the high dimension limit. The high dimension expansion is indeed in 
powers of 1/D and it coincides at the first order in D with the Xjz expansion. 

6 Appendix. 

In order to demonstrate that (|l|) is the correct functional for the free energy we can show that 
lim^^o is the internal energy when is solution of (j^) and that (|l|) is correctly 
normalized at /5 = {—(3f{f3 = 0) = In 2); the latter condition is easy to verify considering 
that = 0) = 1 and that Tr^-^ gives 2"" terms. To convince oneself of the validity of the 
former assertion one can construct explicitly the order parameter gniWo}) making clear its 
physical meaning. Following the approach of we start writing the partition function in a 
recursive manner making use of the equivalence of the model with a Cayley tree. Focusing on 
an arbitrary spin ctq: 

z+l 

Z = J2exp {(3hao) U QiL){(ToW^'^) , (40) 
M fe=i 

where 

Q^L){(ToW^'^) = exp{(3Joiao(T, + (3hai) fl g(L-i)(cr, . (41) 

k=l 

^This is true for the Gaussian model, but it is not true for the ±1 at fixed z, where spin fou are present. 
However it is reasonable that this difference can be seen only at higher orders in the 1 /z expansion. 
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z + 1 is the branches number (random lattice's connectivity) and L the shells number; the a'^^^ 
are the spins on fc-th branch excluding the ctq; h is an external uniform field. 
We can than write the n replicas partition function: 

= n = E ■ ■ ■ E exp E n n Qi^Ma^'^"") m 

a=l {""n} \a=l / fc=l a=l 

and define: 



^7n,(L)(K})= E n%)Kk^'=^") ' (43) 

{^{fc)"}a=l 

where the bar is the average over the random couplings J . 
From and (^) we obtain: 

^= Ee^pfE/5K)^7:ji)(K}) , (44) 

K} ^"=1 ^ 

that reveals 5'n,{L)({o"o}) ^e the one branch contribution to the partition function. 
By definition it follows the recursion relation: 

P + OO / " " \ 

9nXL){{<})= / rfJP(J)Eexp E/^^^ + E/^-^^oX <(L-i)(K})- (45) 

{<7j} Va=l a=l J 



The internal energy density can be written as a bond energy multiplied by the number of 
links per spin {{z + l)/2). 

If we consider a link with a coupling constant J between two spins do and di, we can write 
its energy as 0: 

^01 = < f^o^i > • (46) 

The expectation value is computed with the Hamiltonian H = —Ja^ai + + Hi^ where 
Hq is the Hamiltonian of the spin (Tq before being connected with ai and can be written as 
Hq = — In. [gniWi)}))^ the same argument can be repeated for cxi. 
At this level we should use the finite normalized order parameter: 



which follows the recursion equation: 

/+^^^JP(J)EK}e^P(Sa=l/?^^0"^?)<(L-l)({<}) ,,,, 
t/n,(L)(|(Xo|) = ry— yr . (48) 

In the thermodynamic limit, taking into account only the inner part of the Cayley tree, the gn 
are shells independent, so we can write: 

gniXOa)) = 77]— T\ ' (4^) 
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Anyway the internal energy (as the free energy) is insensitive to the normahzation of the 
order parameter, so we can use both (^9]) or (^51) (in the last one the thermodynamic limit has 
to be taken). 

We can now evaluate the derivative of (3f with respect to /3 and appurate that we obtain 
the expression we expected: 

/_+^ rfJP(J)Tr{.4rr{.4^:({a4)^-({r4) exp [/3JELi ^ar,](-JELi ^g^g) 
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